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We report a first-principles study of a class of (BiScOs^-^PbTiOa)^ (BS-PT) alloys recently 
proposed by Eitel et al. as promising materials for piezoelectric actuator applications. We show that 
(i) BS-PT displays very large structural distortions and polarizations at the morphotropic phase 
boundary (MPB) (we obtain a c/a of ~1. 05-1.08 and P tct ~ 0.9 C/m 2 ); (ii) the ferroelectric and 
piezoelectric properties of BS-PT are dominated by the onset of hybridization between Bi/Pb-6p 
and 0-2p orbitals, a mechanism that is enhanced upon substitution of Pb by Bi; and (iii) the 
piezoelectric responses of BS-PT and Pb(Zri_ a; Ti a; )03 (PZT) at the MPB are comparable, at least 
as far as the computed values of the piezoelectric coefficient dis are concerned. While our results are 
generally consistent with experiment, they also suggest that certain intrinsic properties of BS-PT 
may be even better than has been indicated by experiments to date. We also discuss results for PZT 
that demonstrate the prominent role played by Pb displacements in its piezoelectric properties. 



I. INTRODUCTION 

Perovskite alloys based on PbTiO,3 (PT) are of con- 
siderable interest for applications as piezoelectric actua- 
tor materials. 1,2 The phase diagrams of the most tech- 
nologically important alloys are characterized by a mor- 
photropic phase boundary (MPB) separating the PT-rich 
tetragonal phase, in which the polarization lies along a 
(001) direction, from the PT-poor rhombohcdral phase, 
in which the polarization is along a (111) direction. 3 Be- 
cause the structural transition between these two phases 
brings about a large electromechanical response, materi- 
als with a composition lying close to the MPB are the 
preferred ones for applications. Examples of such mate- 
rials are PZT or PZN-PT, in which PT is alloyed with 
PbZr0 3 (PZ) or Pb(Zn 1/3 Nb 2/3 )03, respectively. 

Recently, attention has been drawn to a new class 
of materials in which PT is alloyed with Bi-based per- 
ovskites, the best-studied example being (BiSc03)i_ x - 
(PbTiOa^ (BS-PT). 4 The appeal of BS-PT is twofold: 
(i) it appears to have piezoelectric properties compara- 
ble to those of PZT and PZN-PT in quality, and (ii) 
its dielectric and piezoelectric properties should be more 
robust, over a wider temperature range, than those of 
PZT and PZN-PT. The physical reason behind (ii) is 
that the Curie temperature within the MPB composition 
range (T^ PB ) is higher in BS-PT (~ 450°C) than in PZT 
(~ 400°C) or PZN-PT (~ 200°C), suggesting that BS-PT 
at room temperature will age more slowly, its properties 
will be more temperature-independent, etc. 2,4 A possible 
drawback is that BS-PT is evidently not thermodynam- 
ically stable in the perovskite crystal structure over the 
entire range of composition. However, for applications 
purposes it is only important that the perovskite struc- 
ture can be obtained for compositions near the MPB, and 
this is indeed the case. 

The authors of Ref. 4 synthesized BS-PT guided by 
an empirical crystal-chemistry rule suggesting that since 



Bi 3+ is too small an ion to form a cubic perovskite with 
Sc 3+ , it will have a strong tendency to move off-center 
from its high-symmetry position, yielding a high struc- 
tural transition temperature. While perhaps overly sim- 
plistic, this reasoning is undoubtedly partially correct 
and is likely to be useful for materials engineering. How- 
ever, it is unclear how to explain the large piezoelectric 
responses found in BS-PT based on ion-size considera- 
tions alone. The ferroelectric properties of materials like 
PZT are known to be related to the partial covalency 
of some bonds, 5 ' 6 and the situation in BS-PT should be 
similar. On the other hand, most of the relevant piezo- 
electric alloys have only Pb on the A site of their per- 
ovskite structure, and the chemistry of Pb is known to 
be very important to their properties. 7 The substitution 
of Pb by Bi in BS-PT is thus of particular interest in 
view of the good piezoelectric properties of this material. 

We present here a first-principles study of BS-PT. We 
find that Bi plays a crucial role, and in particular that hy- 
bridization between Bi-6p and 0-2p orbitals is the driving 
force for the ferroelectric instabilities of BS-PT, allowing 
for very large polarizations and responses. Our results 
are generally consistent with experiment, but they also 
suggest that the intrinsic ferroelectric and piezoelectric 
properties of BS-PT alloys might be even better than 
those measured experimentally to date. 

We present our results in the form of a systematic 
comparison of the properties of BS-PT with those of 
PZT. Section II describes the first-principles methods 
employed. In Sec. Ill we study the ferroelectric insta- 
bilities of pure BS, PT, and PZ, introducing some cor- 
responding calculations on B1YO3 (BY) as an aid to the 
discussion. In Sec. IV we discuss the ferroelectric insta- 
bilities of the alloys and make contact with the experi- 
mental results, while Sec. V is devoted to the piezoelec- 
tric properties of the alloys near the MPB. We summarize 
and conclude in Sec. VI. 
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TABLE I. Relaxed structure and polarization for BiSc0 3 (BS), PbTi0 3 (PT), PbZr0 3 (PZ), BiYO s (BY), a VCA BS-PT 
alloy at x = 0.5 (BS-PT (VCA)), a 10-atom BS-PT supercell (BS-PT (SC)), a VCA PZT alloy at x = 0.5 (PZT (VCA)), and 
a 10-atom PZT supercell (PZT (SC)). (See text for details.) Reported are the tolerance factor t; the energy difference AE (in 
eV) per 5-atom cell relative to the cubic phase, and the polarization P (in C/m 2 ), for tetragonal and rhombohedral phases; 
and the c/a ratio and the rhombohedral angle a (in degrees) for the tetragonal and rhombohedral phases respectively. BS-PT 
(SC) has two inequivalent rhombohedral phases; the reported results correspond to the lower-energy phase. 



System 


t 


AE 


Tetragonal phase 
c/a 


P 


Rhombohedral phase 
AE a 


P 


BS 


0.907 


-1.124 


1.285 


0.93 


-1.353 


88.6 


0.47 


PT 


1.027 


-0.060 


1.049 


0.82 


-0.047 


89.6 


0.72 


PZ 


0.970 


-0.214 


1.035 


0.73 


-0.274 


89.5 


0.78 


BY 


0.845 


-2.260 


1.376 


0.95 


-2.896 


88.7 


0.49 


BS-PT (VCA) 




-0.376 


1.114 


0.92 


-0.412 


89.1 


0.60 


BS-PT (SC) 




-0.472 


1.079 


1.06 


-0.565 


89.4 


1.05 


PZT (VCA) 




-0.060 


1.026 


0.69 


-0.061 


89.6 


0.67 


PZT (SC) 




-0.096 


1.032 


0.73 


-0.104 


89.5 


0.72 



II. METHODS 

The calculations were performed within a plane- 
wave implementation of the local-density approximation 
(LDA) to density-functional theory, using ultrasoft pseu- 
dopotentials 8 to represent the ionic cores. The following 
electronic states were included in the calculation: the 2s 
and 2p states of O, the 3s, 3p, 3d, and 4s states of Sc 
and Ti, the 4s, 4p, Ad, and 5s states of Y and Zr, and the 
5d, 6s, and 6p states of Pb and Bi. For the 5-atom unit 
cell calculations, we used a 35 Ry plane-wave cutoff and 
a 6 x 6 x 6 k-point grid for Brillouin zone integrations. 
We checked that these choices yield converged results. 
We also made some calculations with a rocksalt-ordered 
(fee) 10-atom unit cell, for which we used a 30 Ry cutoff 
and a 4 x 4 x 4 k-point grid. Atomic relaxations were 
considered to be converged when the residual forces were 
smaller than 5 x 10~ 5 a.u. Polarizations were computed 
using the Berry-phase expression of King-Smith and Van- 
derbilt 9 and always using a dense-enough k-point mesh 
along the direction of the strings. The character of the 
electron energy bands was estimated by computing the 
partial density associated with the pseudopotential pro- 
jectors of specific atoms (see Eq. (19) of Ref. 8 for guid- 
ance). 

Calculations using the virtual-crystal approximation 
(VCA) were implemented as in Ref. 10. In the VCA one 
works with virtual atoms defined by a pseudopotential 
that is a weighted average of those of real atoms. Such a 
procedure may lead to meaningful estimates of the aver- 
age properties of an alloy provided that the valence or- 
bitals of the real atoms forming the virtual atom are sim- 
ilar in character. This can be expected to be true for the 
pairs Bi-Pb and Sc-Ti in BS-PT. Also, both groups BiSc 
and PbTi have a total nominal ionic charge of +6, al- 
lowing them to be mixed in arbitrary proportions. These 
requirements are also satisfied by PZT, which has previ- 
ously been studied within the VCA. 10 ' 11 It is important 
to note that the VCA is not expected to render an ac- 



curate description of the electronic structure of a truly 
disordered alloy, but rather to provide us with qualita- 
tive energetic and structural trends. Also, the VCA is 
expected to work better for isovalent pairs of atoms (e.g., 
Zr-Ti) than for heterovalent pairs (e.g., Bi-Pb or Sc-Ti), 
which makes BS-PT a particularly challenging system. 
For this reason, we assess the reliability of our VCA re- 
sults by comparing them with some supercell calculations 
and available experimental information. 

III. PURE COMPOUNDS 

As a preliminary step towards understanding BS-PT, 
we first consider pure BS (in the perovskite structure 12 ) 
and compare with corresponding calculations for pure PT 
and PZ. The first three lines of Table I show the struc- 
tural and polarization results obtained by relaxing these 
systems within ferroelectric phases of tetragonal (Pimm) 
or rhombohedral (R3m) symmetry, always maintaining a 
5-atom unit cell. The atomic relaxations corresponding 
to the tetragonal phases of BS and PT are given in Ta- 
ble II. 

For BS, we find the rhombohedral phase to be lower in 
energy than the tetragonal phase, which agrees with the 

TABLE II. Lattice constant c and atomic coordinates (rel- 
ative to A site, in units of c) in the tetragonal phase of sev- 
eral systems (following the notation of Table I). Note that 
the ideal cubic-perovskite coordinates are z(B)=0.5, z(O x , 
O y )=Q.5, and z(Oz)=0.0. 





BS 


PT 


BY 


(VCA) 
BS-PT 


(VCA) 
PZT 


c (a.u.) 


9.14 


7.65 


10.00 


8.07 


7.65 


z(B) ' 


0.573 


0.535 


0.575 


0.571 


0.547 


2(0*, O v ) 


0.729 


0.606 


0.761 


0.682 


0.610 


z{O z ) 


0.177 


0.094 


0.180 


0.145 


0.085 
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experimental fact that BS-rich BS-PT is rhombohedral, 
and we find that the relaxations are very large. This is 
reflected both in the structural data (we obtain c/a — 
1.29 for the tetragonal phase of BS, compared with 1.05 
for PT; see also the atomic displacements in Table II) and 
in the energetics (AE tc t and AEVho of BS are one order 
of magnitude larger than for PT and PZ). These results 
suggest that the cubic perovskite structure is not a very 
natural one for BS, a conclusion which is reinforced by 
experimental indications that pure BS at zero pressure is 
not even thermodynamically stable. 4 

The polarization associated with the tetragonal phase 
of BS is large, 13 while that of the rhombohedral phase is 
relatively small. In both cases, the largest contribution to 
P comes from the displacement of Bi relative to its neigh- 
boring O ions (O x and O y in Table II). More specifically, 
in the cubic phase of BS all effective charges are normal, 
i.e., with values close to the nominal ionic charges, ex- 
cept for those of Bi (6.4) and and O y (—3.2 along 
the polarization direction). Moreover, these are the ions 
with the largest relative displacements associated with 
the ferroelectric instability (see Table II). The diffcr- 
nce in P between the two phases reflects a significantly 
smaller charge transfer when the Bi ion approaches three 
O ions (rhombohedral phase) instead of four (tetragonal 
phase). 

The fundamental role of Bi in the ferroelectric insta- 
bilities of BS is best appreciated in the electronic band- 
structure data shown in Fig. 1. As in PT (panel c), the 
top valence states of the cubic phase of BS (panel a) 
are 0-2p in character. However, the lowest unoccupied 
states of BS are Bi-6p like, at variance with the Ti-3e? 
bands we find in PT. This fact, together with the rela- 
tively small gap of cubic BS, results in ferroelectric in- 
stabilities mainly characterized by a hybridization be- 
tween Bi-6p and 0-2p orbitals, not the one involving d 
bands of the B cation that is typical of most other ferro- 
electric perovskites. 6 ' 14 ' 15 In the tetragonal phase of BS 
(panel b) , the gap opens considerably and the top valence 
bands increase their Bi-6p character. The strong Bi-6p- 
0-2p hybridization is also reflected in the fact that the 
effective charge of Bi in the tetragonal phase recovers a 
normal value (2.9 along the polarization direction). 

As indicated earlier, Eitel and co-workers chose to 
study BS-PT because of the large size mismatch between 
Bi and Sc. Such a mismatch is conveniently quantified 
in terms of the tolerance factor 

V2{R B +R o y 

where Ro, Ra and Rb are the ionic radii of oxygen and 
of the A and B cations respectively. 16 That is, t is the 
ratio between the ideal cubic lattice parameters based on 
A-0 vs. B-0 bonding alone, with i=l corresponding to 
high perovskite stability. As can be seen from Table I, 
the relatively low value of £=0.907 for BS implies that the 
Bi ion should be very unstable in the high-symmetry po- 
sition of the cubic perovskite structure, consistent with 
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FIG. 1. Near-gap electronic band structure along the T-X 
direction for (a) cubic BS, (b) tetragonal BS, and (c) cubic 
PT. Zero of energy is at the valence-band maximum. 

our findings. We also studied BiYO,3 (BY), for which 
£=0.845, and found even stronger ferroelectric instabili- 
ties that are very similar in character to those of BS (see 
Tables I and II). However, P tct and P r ho of BY are almost 
identical to those of BS, not larger. Also, the distance 
between Bi and its bonded O neighbors is roughly the 
same in the ferroelectric phases of BS and BY (4.33 and 
4.35 a.u., respectively, in the tetragonal phase). These 
two facts seem to indicate that the onset of the Bi-6p-0- 
2p hybridization (which also dominates the development 
of polarization in BY) is not very i-dependent, suggest- 
ing that the large polarizations of BS and BY cannot 
be attributed mainly to a size-mismatch effect. Instead, 
they are strongly linked to the particular chemistry of 
Bi, i.e., its strong tendency to bond covalently with the 
surrounding oxygen anions. 



IV. ALLOYS: FERROELECTRIC INSTABILITIES 

We now turn to a study of the BS-PT alloy within 
the virtual crystal approximation (VCA). 10 ' 11 We also 
carry out parallel VCA calculations on the PZT alloy as 
a relevant reference. 

Figure 2 shows the energies of the tetragonal and 
rhombohedral ferroelectric phases of BS-PT and PZT. 
(The cubic phase, which is material- and composition- 
dependent, is chosen as the zero of energy.) Structural 
and polarization results corresponding to those obtained 
for the pure compounds are shown in Tables I and II. 

Associating the MPB composition with the crossing 
of the energy curves for the two ferroelectric phases in 
Fig. 2, we estimate xmpb ~ 0.65 for BS-PT, while the 
experimental value is 0.64. 4 Similarly, for PZT we get 
Xmpb ~ 0.55, to be compared with experimental val- 
ues from about 0.45 to about 0.52 (the range of stability 
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FIG. 2. Calculated VCA energies of the tetragonal and 
rhombohedral phases, relative to the cubic phase, as a func- 
tion of PT content for BS-PT (panel a) and PZT (panel b). 
Note the two different energy scales. 



of the intermediate monoclinic phase of PZT 3 ). Taking 
into account all the approximations involved in our es- 
timation, we must regard the agreement in the case of 
BS-PT as partly fortuitous. Nevertheless, these results 
suggest that our approach is adequate. 

Figure 2 also shows that the energy reductions asso- 
ciated with the ferroelectric instabilities of BS-PT are 
considerably larger than those of PZT. In particular, at 
the MPB we have A£ MPB = -0.234 eV for BS-PT and 
—0.056 eV for PZT. These results are reminiscent of the 
highly unstable cubic phase of BS previously discussed. 

It is tempting to try to connect these numbers with 
the Curie temperatures of BS-PT and PZT at the MPB. 
We know that, in principle, the thermal energy required 
for the system to remain in the high-symmetry paraelec- 
tric phase is of the order of |Ai? MPB |, which suggests 
that T<¥ PB of BS-PT is about 5 times larger than that 
of PZT, i.e., higher than 1000°C. However, we also know 
that in systems in which there is competition between 
different structural instabilities, such an argument can 
grossly overestimate the transition temperature. In PZT, 
for example, it is known that there is a competition be- 
tween the ferroelectric rhombohedral (FE-rho) phase and 
an antifcrrodistortive (AFD) phase involving rotations of 
the oxygen octahedra. 17 Such a competition also exists in 
BS-PT, since we find that this material presents a strong 
AFD instability. Taking PZT as a reference, one could 
try to estimate the effect of this competition in BS-PT in 
the following way. The AFD phase supposes a strain of 
the high-symmetry structure that is opposite to that of 
the FE-rho phase. At for BS-PT the volume per 

5-atom cell of the FE-rho phase at the MPB is 415 a.u., 
while it is 396 a.u. for the AFD phase, resulting in a vol- 
ume ratio of 1.05. This ratio is 1.01 in PZT, which sug- 
gests that the FE-rho- AFD competition will be stronger 
in BS-PT than in PZT. Energetically, though, we find 
the opposite trend. The FE-rho phase is more favorable 



by 0.15 eV in the case of BS-PT, while in PZT FE-rho is 
preferred to AFD by only 0.02 eV. Hence, we are unable 
to determine in which of the two systems the reduction of 
T^ PB caused by the FE-rho-AFD competition should be 
larger. On the other hand, there is another mechanism 
missed by the VCA and which could significantly lower 
the transition temperatures of BS-PT. Because BS-PT is 
a disordered heterovalent alloy, internal electric fields will 
be present in it. Such fields will cause local frustrations 
and enhance competition between different polarization 
directions, resulting in a decrease of T^. 18,19 In isovalent 
PZT the internal fields and their effects are bound to be 
much smaller. Quantifying these effects would require a 
full statistical-mechanics study of the disordered alloys, 
but that is beyond the scope of this work. In conclusion, 
all we can say about the issue of the transition temper- 
atures is that our results indicate that J* IPB of BS-PT 
could well be significantly larger than that of PZT. 

Let us now discuss the structural and polarization 
data. The ionic and strain relaxations of BS-PT arc 
much larger than those of PZT (see Tables I and II). Im- 
portantly, at the MPB composition we find c/a =1.079 
and P tct = 0.875 C/m 2 for BS-PT, to be compared, re- 
spectively, with 1.027 and 0.687 C/m 2 obtained for PZT. 
The very large values calculated for BS-PT in the MPB 
region suggest that its piezoelectric response could be 
very large too, assuming the usual picture in which this 
arises from easy polarization rotation associated with the 
near-degeneracy between tetragonal and rhombohedral 
phases. 20 ' 21 In fact, according to our results, it would 
not be surprising to find that BS-PT has even better 
piezoelectric properties than PZT. 

Here we have to note that, while our results for PZT 
are in reasonable agreement with experiment (at the 
MPB, c/a and P of PZT are measured to be about 1.025 
and 0.75 C/m 2 , respectively 22 ' 23 ), our results for BS-PT 
significantly deviate from experimental values. Ref. 4 re- 
ports c/a w 1.023 and P tct w 0.32 C/m 2 at the MPB 
of BS-PT, values that are much smaller than those that 
result from our VCA calculations. In view of this dis- 
crepancy, we decided to asses the reliability of the VCA 
by studying a 10-atom superccll of BS-PT at 50/50 com- 
position. We have chosen an fee superccll in which there 
is rocksalt ordering both on the A sublattice (Bi/Pb) and 
on the B sublattice (Sc/Ti). The resulting supercell has 
Td symmetry so that all tetragonal ferroelectric domain 
orientations remain degenerate, while the the rhombohe- 
dral domains split in two groups of four degenerate ori- 
entations. (Here we report results for the lower-energy 
rhombohedral domains only.) The results, summarized in 
Table I, display a significant reduction in the c/a of the 
tetragonal phase, but not large enough to reconcile the- 
ory with experiment. The VCA error is presumably max- 
imum at x = 0.5, so at Xmpb we would expect that c/a 
should be at least 1.045, still significantly larger than the 
experimental value. Regarding the polarization, the 10- 
atom supercell values are actually larger than the VCA 
results. 
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In Table I we also included results corresponding to 
the analogous 10-atom supercell calculations for PZT. As 
expected, 10 the VCA is more accurate for PZT (isovalent 
alloy) than for BS-PT (heterovalent) . This is further con- 
firmed by the atomic relaxation data (not shown here). 

Hence, for BS-PT a relatively large discrepancy be- 
tween theory and experiment remains. Its origins arc 
most probably related to the differences between the ide- 
alized crystalline samples we consider and the ceramic 
samples actually grown (the internal fields present in real 
disordered samples are likely to play a particularly rele- 
vant role). In any case, we hope that the very large cja 
and polarizations we obtain at the MPB will stimulate 
further experimental study of this promising system. 24 



TABLE III. Effective charges Z\^ and derivatives dui/dr/s 
associated with the eis,i piezoelectric coefficients of VCA 
BS _p T *=o.7 and pZT x=o.6 Derivatives given in units of a, 

the lattice parameter along the x direction; a— 7.27 a.u. for 
BS _ PT *=o.7 and a=7 42 a . u . f or PZT *=o.6 We uge the CQn _ 

vention that dui(s)/drjz = 0. 



dui /dr/5 





BS _p T *=0-7 




BS-PT^ 07 




A 


4.26 


3.82 


1.00 


0.92 


B 


5.27 


6.64 


0.07 


0.16 


o x 


-4.47 


-5.23 


0.09 


0.07 


Oy 


-2.77 


-2.90 


-0.41 


-0.51 


o z 


-2.29 


-2.33 


-0.76 


-0.64 



V. ALLOYS: PIEZOELECTRIC RESPONSE 
NEAR THE MPB 

Finally, we compute the piezoelectric response of BS- 
PT, and compare with that of PZT. First we investigate 
the coefficients = dPi / dr\j that describe the so-called 
inverse piezoelectric effect (r/j is the j-th component of 
the strain tensor in Voigt notation). 26 We report here 
our results for ei 5 of the tetragonal phase, the coeffi- 
cient known to be responsible for the large piezoelectric 
response of PZT. 21,27 We focus on VCA alloys close to 
the MPB on the tetragonal side; specifically, we consider 
BS-PT at x = 0.7 (BS-PT 2 ^ - 7 ) and PZT at x = 0.6 
(PZT x=oe ). 

Following Ref. 28, we split the piezoelectric coefficient 
into two parts 



eis = ei5 i + ei5 , 



(2) 



where e 15 c is computed at fixed internal atomic coordi- 
nates ("clamped-ion" term) and 



(3) 



("internal-strain" term). Here s runs over the atoms in 
the unit cell, e is the magnitude of the electron charge, 
Z\ x is the effective charge associated with displacements 
u\ along the x direction, a is the lattice constant in the 
x direction, and v is the unit cell volume. 

For BS-PT 2 ^ - 7 wc obtain ei 5 , c = -0.93 C/m 2 
and ei5 ; j = 8.18 C/m 2 , which add up to eis = 
7.25 C/m 2 . For PZT 2 ^ 06 we obtain ei 5 , c = -0.66 C/m 2 , 
ei5 ; i =8.24 C/m 2 , and e 15 = 7.58 C/m 2 . Hence, both sys- 
tems have a similar ei 5 , and in both the internal-strain 
contribution strongly dominates the response. To under- 
stand this, we note that a large value of eis.i can arise 
in Eq. (3) cither from a "dielectric effect" (large Z*) or 
an "elastic effect" (large response of the internal coordi- 
nates to strain, i.e., large dui/dr/^). Table III shows the 
calculated values of Z\ x and du\jdr\^ for all atoms in the 



unit cell of VCA BS-PT X = 7 and PZT X = - 6 . In both sys- 
tems the response is mainly associated with an enhanced 
covalent bonding between the A ion and oxygens O y and 
O z , as pictorially shown in Fig. 3. Our results further 
suggest that the substitution of Pb by Bi enhances such 
a mechanism, consistent with the fact that Bi is more 
electronegative than Pb. The B ion contributes by shift- 
ing closer to O^, but for the B-Oa; pair the elastic effect 
is substantially smaller than it is in the case of A-O y 
and A-O z . We also note that when these results are 
compared with the analogous ones for the smaller e 33 co- 
efficient (not shown here - see Refs. 10 and 28 for related 
PZT results), one finds that it is the elastic effect that 
makes eis particularly large. 

The above results on eis have revealed much about 
the chemical origin of the piezoelectric response of these 
systems. However, the energetics of the piezoelectric re- 
sponse, as well as the actual technological interest of a 
material, is best characterized by the direct piezoelec- 
tric coefficients dij = dPijdoj reflecting the polarization 
that develops in direction i by application of a stress with 
Voigt component j (or cquivalently, the strain j induced 
by an electric field along i). In particular, dis measures 
the ease of polarization rotation associated with the large 
responses at the MPB. Technical difficulties have tradi- 
tionally hampered the direct calculation of dij coefficients 
from first principles. 29 In the present case, however, it 
is not hard to show that di 5 is just given by ei 5 /C 5 , 
where C5 is the shear elastic constant associated with 
775 when the atomic positions are allowed to relax in re- 
sponse to the strain (i.e., the dressed elastic constant) 
at fixed electric field. We calculated g?is =168 pC/N 
(with l7 5 =0.62 a.u.) for BS-PT x=a7 and d 15 =196 pC/N 
(with C 5 =0.56 a.u.) for PZT X=0S . These results agree in 
magnitude with the piezoelectric coefficients of ceramic 
samples of BS-PT and PZT measured at room temper- 
ature and at the MPB, as well as with the semiempiri- 
cal calculation of Wu and Krakauer for an ordered PZT 
supercell. 27 

Thus, our calculations clearly provide evidence that 
BS-PT can be expected to be a piezoelectric of com- 



5 



parable quality to PZT. This conclusion, as well as the 
leading role of Bi and/or Pb in the response, is further 
supported by calculations of e 33 in the tetragonal phase 
of the above-mentioned 10-atom supercells (i.e., by cal- 
culations in which the VCA was not used). 

VI. SUMMARY AND CONCLUSIONS 
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Our first-principles study of BS-PT indicates that this 
material displays very large structural distortions and po- 
larizations at the MPB. In particular, we obtain c/a be- 
tween 1.05 and 1.08 and P tct w 0.9 C/m 2 . We also find 
that the piezoelectric response of BS-PT near the MPB 
is comparable to that of PZT. 

Our calculations show that the large polarization and 
piezoelectric responses in BS-PT are mainly related with 
the onset of Bi/Pb-6p-0-2p hybridization, a mechanism 
that is enhanced upon substitution of Pb by Bi, since 
Bi is a more covalent atom. Our results also provide 
evidence of the dominant role of Pb in the piezoelectric 
response of PZT, in agreement with recent experimental 
studies. 7 In both BS-PT and PZT, this predominance of 
the A-site ions is related to a large elastic effect, i.e., a 
large displacement of the ions in response to strain. 

Our results are in reasonably good agreement with ex- 
periment. The quantitative discrepancies between theory 
and experiment (in particular, for the values of c/a and 
the polarizations at the MPB) suggest that the intrinsic 
piezoelectric properties of BS-PT alloys may be even bet- 
ter than those measured experimentally to date. We thus 
hope that our results will stimulate further experimental 
work on this and related materials systems. 
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FIG. 3. Sketch of the unit cell of the tetragonal phase of 
VCA BS-PT or PZT. The arrows represent the displacements 
associated with eis = dP\/dr\h- Atoms in the pairs A-O y , 
A-O z , and B-Oa, approach each other. 
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